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Abstract 

We mainly study the Renyi entropy and entanglement entropy of the states locally 
excited by the descendent operators in two dimensional conformal field theories (CFTs). 
In rational CFTs, we prove that the increase of entanglement entropy and Renyi entropy 
for a class of descendent operators, which are generated by onto the primary 

operator, always coincide with the logarithmic of quantum dimension of the corresponding 
primary operator. That means the Renyi entropy and entanglement entropy for these de¬ 
scendent operators are the same as the ones of their corresponding primary operator. For 
2D rational CFTs with a boundary, we confirm that the Renyi entropy always coincides 
with the logarithmic of quantum dimension of the primary operator during some periods of 
the evolution. Furthermore, we consider more general descendent operators generated by 
(Hi L-rti rij L-rij ) on the primary operator. For these operators, the entangle¬ 
ment entropy and Renyi entropy get additional corrections, as the mixing of holomorphic 
and anti-holomorphic Virasoro generators enhance the entanglement. Finally, we employ 
perturbative CFT techniques to evaluate the Renyi entropy of the excited operators in de¬ 
formed CFT. The Renyi and entanglement entropies are increased, and get contributions 
not only from local excited operators but also from global deformation of the theory. 
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1 Introduction 


The entanglement entropy (EE) and the entanglement Renyi entropy (RE) are very helpful 
quantities to study global or non-local structures in quantum field theories. For example, 
computing topological entanglement entropy can characterize the topological order of the 
system |lj. In [2] a relationship between the topological entanglement entropy and boundary 
entropy has been pointed out. And in [3j the relationship between the boundary entropy and 
entanglement entropy has been explored. Both the entanglement and Renyi entropes of a 
subsystem A are defined with respect to the reduced density matrix pa, which is obtained by 
tracing out the degrees of freedom in the complement of A from the original density matrix 
p. Formally, the n-th Renyi entanglement entropy is defined by S^ = log Tr[//|]/(1 ~ n). 
If the limit n —>• 1 is well-defined, lirn S\ ' coincides with the von-Neumann entropy exactly, 

n—>i 

which defines the entanglement entropy. This leads to standard replica trick to in compute 
the entanglement entropy in the quantum held theory. 

There could be topological contribution to the entanglement entropy even for gapless the¬ 
ories, e.g. conformal held theories (CFTs). In 2D rational CFTs, it was found [II] that for the 
locally primary excited states, the Renyi entropy difference is related to the quantum dimen¬ 
sion a a of the primary operator, which is a kind of topological quantity. The computations 
of the entanglement entropies for locally excited states have been formulated in mmm in 
the field theory side. The entanglement entropy for local free scalar have been investigated 
in m 0 (see also [10] for the operator approach). In the large N CFTs , the entanglement 
entropies for locally excited states have been discussed in [MEO- In a more recent paper m, 
the Renyi entropy difference Aof the local excited states in a boundary conformal held 
theory (BCFT) has been studied. In [16] the left-right entanglement entropy for a boundary 
state has been discussed. 

In all the study of quantum entanglement of local operator, only the primary operators 
have been discussed^ It is interesting to investigate the similar effect from a local descendent 
operators. Among the descendants, the stress tensor is of particular interest, as it is related 
to the energy density in the held theory, and is dual to a graviton in the holographic dual. 
In this paper, we generalize the previous study ffn m on the Renyi entropy difference to 
the descendent operators with or without a boundary condition. First of all, we compute 
A<!v 2 ' explicitly for some simple descendent operators in minimal models, and we find that 
in these cases the entropy difference is exactly the logarithmic of the quantum dimension of 
corresponding primary operator. Moreover, we study the descendant operators in a general 
theory. The picture is that if the local operator is of the form the entropy difference 

1 There has been some discussion on the entanglement entropy for descendant excited states in ]7]. 
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is the same as the one for the local primary operator. This could be explained by studying 
the divergent terms in the correlation function. However for a more generic operator like 

V= d {n i },{n j }Yl L -n i l- nj O(w,W ) ( 1 . 1 ) 

{ni},{nj} i,j 

we find that the Renyi entropies and entanglement entropy have additional contribution. Such 
additional correction originates from the entanglement in the excited states. We give a sys¬ 
tematic recipe to calculate the entanglement entropy for general descendant operators. 

We also investigate the Renyi entropy for the states excited by a local descendent operator 
in the rational CFTs with a boundary. The boundaries introduced here do not break the 
conformal symmetry. Such kinds of theories are called the boundary conformal field theories 
(BCFT). We find that the Renyi entropies excited by the descendants in BCFTs are similar 
to the one discussed in mi jl5j . If the descendent states are generated by the operator (1.1), 
the boundary does change the time evolution of the Renyi entropy but it does not change the 
maximal value of the Renyi entropy. Finally, we investigate the Renyi entropies in a deformed 
CFT perturbatively. If the deformation can be seen as a perturbation, the Renyi entropies of 
locally excited states are just the summation of contributions from the local excitation and 
the one from global deformation. 

The layout of this paper is as follows. In section 2, we would like to study the Renyi 
entropy of the local descendent state with or without a boundary. In section 3, we move to 
study the Renyi entropy of the local excited states in 2D CFT with deformation by chemical 
potentials. In these deformed theories, we obtain the Renyi entropy of a subsystem with time 
evolution. In section 4, we devote to the conclusions and discussions. In appendix, we list 
some technics which are very useful in our analysis. 

While proceeding with this project, we noticed a recent papers [3J] appearing in arXiv, 
which has overlap with our discussion in the sections 2.3 and 2.4. 

2 Renyi entropy of descendent operator 

2.1 Setup in 2D CFT 

Consider an excited state which is defined by acting a primary or a descendent operator V a 
on the vacuum |0) in a two dimensional CFT. We make use of the Euclidean formulation and 
introduce the complex coordinate (w, w) = (x + it. x — it) on R 2 such that r and x are the 
Euclidean time and the space respectively. We insert the operator V a at x = — l < 0 and 
consider its real time-evolution from time 0 to t under the Hamiltonian H. The corresponding 
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density matrix can be expressed as following: 


pit) = Af-e- iHt e- eH V a i-l)\0)(0\V^{-l)e- eH e iHt 

= M ■ V a (w 2 , u7 2 )|0)(0|V r a t (i(;i,u)i), (2.1) 

where M is fixed by requiring Tr p(t) = 1, and 

w\ = i(e — it) — l, u >2 = — i(e + it) — l , (2-2) 

wi = —i(e — it) — l, W 2 = i(e + it) — l. (2-3) 


The infinitesimal positive parameter e is an ultraviolet regularization factor. We treat e ± it 
as purely imaginary numbers until the end of the calculations, as in [8l(9l!l7j. 

To calculate A S\, we can apply the replica method in the path-integral formalism by 
generalizing the formulation for the ground state [3] to our excited states 0. It leads to a 
n-sheeted Riemann surface E n with 2 n operators V a inserted. We choose the subsystem A to 
be an interval 0 < x < L at r = 0. In section 2.2, 2.3, we would like to take L 
up. For later parts, we will consider the situation with finite L. 

Finally, the Acan be computed as 

: (n) _ 1 


oo as a warm 


A5\ 


1 — n 


log{V^(wi,W 1 )V a (w2,W 2 ) • • • V a {w2n,W2n) 


-nlog(v£(wi,wi)V a {w2,w 2 ) 


Si 


(2.4) 


where (w 2 k+i, W 2 k+ 2 ) f° r k = 1 , 2 ,..., n — 1 are n — 1 replicas of (wi,w 2 ) in the k -th sheet of 
T, n . The term in the first line is given by a 2n-point correlation function on £ n , while the 
one in the second line is given by n two-point functions on Si. Here A a is the (chiral and 
anti-chiral) conformal dimension of the operator O a . 


2.2 Convention 

Let us study the second Renyi entanglement entropy in details to set up the convention. The 
( 2 ) 

calculations of Ais related to the four-point functions which we know pretty well for 
exactly solvable CFTs. 

A (2) 

To compute A S y A we need to introduce the coordinate Z{. which is related to Wi by the 
conformal map wi = z'f. First the coordinates Wi and Zi behave like 

wi = ie + t — l = re 101 = (^i) 2 , 
w 2 = -ie + t — l = se 102 = {z 2 ) 2 , 

W3 = ( ie + t- l)e 27ti = re^ 2n+0l) = {z 3 ) 2 , 

W4 = {-ie + t - l)e 2ni = se^ 27r+02) = (z 4 ) 2 . (2.5) 
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Figure 1: This figure is to show our setup in two dimensional plane u = x + it . The system 
will be triggered at x = — l and there are left- and right-moving quasi-particles at t = 0. 


Thus we find 


z\ = — Z 3 = y/wi = \/re ldl / 2 = i\Jl — t — ie, 

Z 2 = —Z 4 = y/w 2 = \fse q ’ 9 ' 2 / 2 = iVl — t + ie. (2-6) 

The most important point is that we should treat ±ie +1 as a pure imaginary number in all 

algebraic calculations and take t to be real only in the final expression of the entropy. This is a 

standard prescription to make analytical continuation of Euclidean theory into its Lorentzian 
version. Therefore we should identify 

(r cos 0i, r sin 9\) = (— l, e — it), (s cos 62 , s sin 62 ) = (— l, — e — it). (2-7) 

This leads to 


r = \Jl 2 + (e - it) 2 , s = \Jl 2 + (-e - 
rs = \J ( l 2 + e 2 — t 2 ) 2 + 4e 2 t 2 , r 2 + s 2 
cos($i — O 2 ) = 2 cos 2 1 = 


it) 2 , 

= 2{l 2 + e 2 — f 2 ), 
l 2 - e 2 - t 2 

y/(l 2 + e 2 - t 2 ) 2 + 4 e 2 t 2 ' 


( 2 . 8 ) 


To evaluate the four point-correlation functions, it is useful to focus on the ratio 


^ 12^34 ~(l-t) + y / (^-f) 2 + e 2 

z —-=- —— - 

^13^24 2 y/(l - t ) 2 + e 2 

^ 12^34 —(l + t) + \J(l + t) 2 + e 2 

z -=- - 

^13^24 2^/(1 +t) 2 + e 2 

where Zij = Zi — Zj. Note that there is a useful relation 


^ 14^23 
213^24 ’ 


( 2 . 10 ) 
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We are interested in the two limits (i) l » t » e (early time) and (ii) t » l » e (late 
time). From (12.91) we find that they correspond to 

e 2 

(0 z - z -y2^ °)> 

(**)*- 1 “ ^2 (-*!)’ ^-^2 (^°)' ( 2 ' n ) 

Note that the second limit is non-trivial in that it does not respect the complex conjugate 
relation because of our analytical continuation of t. 


2.3 Second Renyi entanglement entropy in terms of 4-pt functions 

Let us first consider the Renyi entropy of the simplest descendent operator L_iO, where O is 
a primary field. The four-point function on n-sheeted Riemann surface E n can be mapped to 
the one on R 2 by the conformal map w = z n . For the second Renyi entropy, we find 


{L-i 0 (w\,w\)L-\ 0 (w 2 , w 2 )L~iO(w 3 , w 3 )L_ 1 0 (w 4 , w 4 )} 


s 2 


n 


i =1 


dwj 


dzi 


( -A 


dz\ d 2 w\ dz\ 


dwi dz? dw\ 


+ 0(z u z{) -A 


dz 2 d 2 w 2 dz 2 


dw 2 dz% dw 2 


-A 


dz 3 d 2 w 3 dz 3 


dw 3 dz\ dw 3 


+ 0(z 3 ,z 3 ) -A 


dz4 d 2 W4 dz4 


dw4 dz% dw4 


+ 77— d Zi 0 (z 4,Z4 ))t : 1 


= 2 8A \ziz 2 z 3 Z4\ 2/1 • A1A2A3N4I A 4 (0 (zi,z 1 )0 (z 2 ,z 2 )0 (z 3 ,z 3 )0 (z 4,Z4 ))i 


I -2 A 




M, 


A 3 ^2 j2 d zi{0(zi,zi)0(z 2 , z 2 )0(z 3 , z 3 )0(z 4 , Z 4 ))z 1 


i=l 

4 


a 2 x—^ MiMj 

+A 2 2^ 3 


N i N j “* ^ 

^ / j 


d Zi d z . (0(^i, zi)0(z 2 , z 2 )0(z 3 , 23)0(24, 2 4 )) Si 


' A M 2^T^T k d Zi d Zjdz h ( 0 (Z!, Z 1 ) 0 (z 2 , Z 2 ) 0 (z 3 , 23)0(24, 24)) Sl + 


i^j^k 
Ml M 2 -M 3 Ad 4 


NiNjNk 

d zl d Z2 d Z3 d Z4 (0(zi , 21)0(22,22)0(23,23)0(24,2 4 )) Sl 


A1A2A3A4 

where A is the chiral conformal dimension of the primary operator O, and 


A,; = 


dz, d 2 


Wi 


dvjj dz 2 


Mi = 


dzj 

dwi 


+ 0(22,22) 


( 2 . 12 ) 


(2.13) 


Since A and M, t do not contribute to the divergent terms in both the late time limit and the 
early time limit, we can take them as a factor to simplify our analysis. In (12.121) . the last term 
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in the summation diverge fastest in the small e limit, so we only keep this term in the later 
calculation. Note also that rs behaves in the two limits as 


(z)Early time : rs ~ l 2 , 

(ii )Late time : rs ~ t 2 . 

Due to the conformal symmetry, the four-point function on R 2 can be expressed as 

(0(zi,Zi)0(z2,Z2)0(z3,Z3)0(Zi,Zi)}^ 1 = |213^24| _4A ' G(z,Z ), 

where (z, z) are given by m- 
The two-point function is 


(2.14) 


(2.15) 


{L-iO(wi,wi)L-iO(w 2 , w 2 )} Si — d wi d u 


Af 2A(2A + 1)A f 


(2.16) 


2 ki 2 | 4A (2e) 4A + 2 

where Af is the normalization of the primary operator O. Note that the four-point function is 
proportional to Af 2 and the Ais of course independent of Af. 

In the late time limit (ii), we finally find that the ratio Z 2 /(Z\) 2 in (12.41) is expressed in 
terms of the four-point function on R 2 : 

2 _ (L- 1 0(wi,Wi)L- 1 0(w2,W2)L-l0{w3,W 3 )L-i0(w4,W4 : )}x 2 


Tr p\ = 


({L_iO(w 1 ,wi)L_iO(w2,w 2 ))z i y 

l 6 e sA + 4 4 


(2A) 2 (2A + 1)W 2 ' —^ N i) d zi^2 d z 3 9z4( 0 ( z i, zi)0(z 2 ,Z2)0{ z 3, z 3 )0(z^ za))^ 

(2.17) 


(2A 2 )(2A + 1)W 2 (4f 2 ) 4A 


1C 8A+4 4 

(II Ni)d Zl d Z2 d Z3 d Zi G(z, 

i =1 


2.3.1 Example I: (2, 1) operator in minimal model 

Now we would like to study an explicit examples: the (2,1) operator in the minimal models. 
Let us consider a (p,p') minimal model with p > p ', where the primary fields 4>(r 1} r 2 ) are 
specified by a pair of integers (ri,r 2 ) taking the values 

n = 1,2,...,p' - 1, r 2 = 1, 2, ...,p — 1. (2.18) 


The central charge of (p,p') minimal model is 

,{p-p’) 2 


c = 1 - 6- 


pp 


j > 


and the conformal dimension of the {r\,r 2 ) primary operator is given by 

(pri — p'r 2 ) 2 - (p-p 1 ) 2 


^ri,r2 — 


4 pp' 


(2.19) 


( 2 . 20 ) 
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Let us focus on the (2,1) operator 4 > (2 ,i) m a (p>p') minimal model as it has a relatively 
simple four-point function. It has conformal dimension 

A «.= fH - (2 ' 21) 

The function G(z,z) for the four-point function (02,1(21, zi)...fa, 1 ( 24 , ^4 ))ei via the relation 
(12.1511 is known to be [THj [T9] 


G(z,z) 


= IzIp 7 !! — z^ 


^ .» SU. (^) i(a 3 + ^ (?) -n (?) a(g 


si ” M 


sin (^) 


where we have defined m m 


( 2 . 22 ) 


h(z) = —^^ • F 

r f 


P 1 . 3p 2p 

~7) — 1 H—7i —; z 

p p' p' 


1 _2P 

J 2 (*) = * p ' 


r 1-5 r 1-5 


r(2-f 


• F 


P , P „ 2p^ 

z: / ■ 2 «> ■2- 

p p' p 


(2.23) 


with F[a,b,c,z] being the hypergeometric function. By taking the limit z \2 = £34 —> 0, we 
hnd that 


(02.1 ( z h Zl) ■ ..02,1 (Z4i ■ 24 ))Ei ->• |^121 


- 8i !nfl( r (‘ ~ 7 ) 1 

“■‘( t e ) r ( 2-^) 2 


(2.24) 


Thus we can identify the normalization factor Af of the two-point function as follows 


M 2 = 


sin(^) 2 rfi-4 4 


sm 


{¥) rb-fF 


(2.25) 


Now we turn to the most interesting limit: late time limit t » l » e, where we have 

2 2 _ 

z = 1 — and z = ^W. Under this limit we can put (12.2211 into the final formula (12.1711 . 

Notice that the I 2 term is more divergent in the small e limit, and using the identities on the 

hypergeometric functions, we have 


h{z) — (1 - z) p 




2 V 


2 \ ^~ 2 vtv' 


Hz, 


~ Z P - 


r( 2 -$ 


FI? 


^ r ( 1_ f) 2 _ / e 2 


r(i 


\4i 2 


2 \ 1 - 2 p/p' 


T 2 


_ p_\ \4t 2 


(2.26) 
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Thus we find 


sm 


G(z, z) 


\: 
ZE P 1 
P' ) 


r (i — ^) r (—1 + |e 


{¥) 


In this way, we finally obtain 


^Pa = 


r(2-| 


r[2-f 


r(-i + f 


2 \ 2-3p/p' 


e 

4t* 


(2.27) 


r (i - 7 ) r 


2 cos 


(2.28) 


and the Renyi entanglement entropy difference 


AS? = log 


o I n P 
—2 cos — 


(2.29) 


For a general minimal model, the entropy difference is not vanishing for the descendent operator 
L_iO(2,1) operator. But for the Ising model ( p,p') = (4,3) [201 [2T] , it is vanishing. 

In a minimal model, the modular transformation maps the primary operator (si,S2) into 
(ri,r 2 ). Its S-matrix is defined to be ‘S , (r-i,r 2 ),(si,s 2 b gi yen explicitly by 


& 

l ^(ri,r 2 ),(s 1 ,s 2 ) = ^ 


l+r 2 si+ris 2 


Sill 


pp 


7Tp 

—ns, ) sm 


irp 


-T 2 S 2 


The quantum dimension d( rijr2 ) for the primary field (ri,r 2 ) is defined by 


i(ri,r 2 ) 


’(i.iMri.rz) 


5, 


o I wp 
= —2 cos — 


(hi), (hi) 


P 


(2.30) 


(2.31) 


Thus we can conclude that for the L_iO(2,1) operator 


A S? = log d (riir2) , 


(2.32) 


the same as the one for the primary operator 0(2,1). 


2.3.2 Example II: energy momentum tensor 

Another simple case is the excitation of energy momentum tenor, which is the descendent state 
of the identity operator. Different from the case discussed above, it is L _2 rather than L_i. 
We only consider AS 2 here. The two-, three- and four-point correlation functions on the R 2 
depend only on the central charge c [22], which are respectively 

(T( Zl )T(z 2 )) = -V, (2.33) 

ZZ 12 


(T( Zl )T(z 2 )T(z 3 )T(z 4 )) 


^12^34 


(2.34) 
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with 


J(z) "i ( 1 + z 4 + (r^) +2c 


z 2 (l-z + z 2 ) 


(2.35) 


(1 -z) A ) (1 -zf 

z = 212234 /( 213224 ) is the ratio. The transformation of the energy momentum tensor under 
the map z(u>) is given by 

dz \2 


rM = (£) T (z) + -^M, 


(2.36) 


with the Schwarzian derivative 


(2.37) 

In the coordinate (12.51) the leading order of the two-point function is 

(TMTM) = (2.38) 

After the transformation the four-point correlation function in the u coordinate is 

r j 7 . 

(T( Ul )T(u 2 )T(u 3 )T(u4)) = H(^-) 2 (T( Zl )T(z 2 )T(z 3 )T(z 4 )) + ..., (2.39) 

i 

where ‘...’ denotes the less divergent terms, such as the product of three-point function and 
the Schwarzian derivative. During the time 0 < t < l or t > L + l, we can see z \ 2 ~ 234 ~ O(e), 
also 2 ~ 0(e 2 ). One could see from (12.341) the leading contribution of the four-point function 
is 

(T( Zl )T(z 2 )T(z 3 )T(z 4 )) ~ (2.40) 

4x -12^34 

with 


iLe 


iLe 


z 12 = -- 


{t - L)(t - L - L) Zu Z34 (t-L(t-L-l ))' 

During the time l<t<L + l, z 1 = — z 2 , 23 = — 24 , 223 ~ 214 ~ 0(e), and 2 ~ 1 — 0(e 2 ), with 

L 2 e 2 


(2.41) 


2 ~ 1 — 


(2.42) 


4 {t-l) 2 (t-L-l) 2 ' 

Now the leading order contribution of the four-point correlation function in the limit e —>• 0 is 

c 2 1 


(T(Z 1 )T( 2 2 )T( 2 3 )T( 2 4 )> ~ 


4(1 - 2) 4 (4212 3 ) 4 ' 


(2.43) 


One could check (12.401) is the same as (12.431) . Therefore, by using (|2.4|) . we find AS^ 2 * = 0 
during the evolution. This is just as we expected, since T(z) is the descendent operator of 
identity I, the quantum dimension of which is 1 . 
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2.4 2nd Renyi entropy for some descendant operators 


In previous section, we computed the Renyi entropy of some specific descendent operators and 
found that these descendant operators had the same contribution as the primary operators. 
In this subsection, we use the conformal block and operator product expansion (OPE) to show 
this is true for a large class of descendent operators. It is also a warm up for the study of more 
general descendant operators in the next subsection. 

In this and next subsection, we will consider the case that the interval is finite, lying at 
[0, L \, and the operators is inserted at —Z. The two-point and 2n-point functions on single 
sheet and n-sheeted Riemann surfaces are respectively 


Z\ = (V(w,w)V(w ,w )), 

Z n = (V(w ll W 1 )V(w[,w[)V(w 2 ,W 2 )V(w' 2 ,W 2 )...V(w n ,W n )V(w n ,W n )) \n-sheets, 


(2.44) 

(2.45) 


where the correlation function Z n is defined on the n-sheeted surface, with V(wi,Wi ) on the 
i-tfi surface. We set 


Wi = —l + (t — ie), Wi = —l + (t + ie), (2.46) 

= — Z + (t + ie), w i = —Z — (t T ie). (2.47) 

V(w,w) can be any descendant operators]^ We only need to set n = 2 for S ^ calculation. To 
evaluate the multi-point correlation function on n-sheeted surface, we need to take a conformal 
transformation 


z = 


w 


y w-L’ 


(2.48) 


First we consider the descendent operator V(w,w) = dO a (w,w), where O a (w,w) is a 
primary operator, the two-point function equals to 


Z\ = (dO a (w, w)dO a (w',w')) 
and the four-point function is 


—2h(2h + 1) 


1 1 

(w — w') 2h+2 (w — w) 2h 


2 h(2h + 1) 
(2e) Ah+2 ’ 


(2.49) 


Z2 = d Wl d w > d W 2 d w ' (0 a (wi,wi)0 a (w 1 ,w 1 )0 a (w2,w2)0 a (w 2 ,w 2 )). 
Under the conformal transformation (12.48(1 . for t < l 

(l — t + ie) 2 _ (l + t — ie) 2 


zi = 


zi = 


{l -\- L — t -{- ie) 2 

(Z — t — ie) 2 
(Z + L — t — ie) 2 


zi = 


*1 = 


(T T l T t — ie) 2 

(Z + t + ie) 2 
{l L t ie) 2 


(2.50) 

(2.51) 

(2.52) 


2 For convenience the convention in this section and the next section is a little different from the one in other 


sections. The correspondence is W2j —> Wj, 102,-1 


to,- 
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Z 2 = 


*2 = 


(l — t + ie) 2 _ (l + t — ie) 2 

— ^2 r 

(/ “h Z/ — t -h ic) 2 (Z/ t — 2£) 2 


(l — t — ie) 2 


^2 = 


(Z + t + ie) 2 


(i T Z/ — t — ie) 2 (/ -f- L -)- t T ie) 2 

where ( 21 , z ^) ( 2 : 2 , -4) (^ 1 , 4) (^ 2 , 4) are close to each other. While for l < t < l + L, 

7ri (t — l — ie) 2 
Zi = -z 2 = e 2 


(2.53) 

(2.54) 


[l -\- L — t T ie) 2 


/ / _e± (t — l 4* ie) 2 

*i = -Z2 = e 2 


(l + L — t — ie) 2 

then (zi,z 2 ) (£ 2 , 4 ) ( 21 ,%) (^ 2 , 22 ) are close to each other. 
From the result in nu, for t < l, 


(2.55) 

(2.56) 


+ less divergent term, (2.57) 


+ less divergent term. (2.58) 


(' O a (zi , zi)O a (z 1 , z 1 )O a (z 2 , z 2 )O a (z 2 ,z 2 )) 

1 1 
(Z\ - z[) 2h (zi — Z^yb (z 2 - z 2 ) 2 h (z 2 ~ z 2 ) 2h 

and for l < t < l + L, 

{O a {zi,zi)O a (z' 1 ,z 1 )O a (z 2 ,z 2 )O a (z 2 ,z' 2 )) 

1 1 1 

da (zi — z 2 ) 2 h {z\ - zy 2h (z 2 - z[) 2 h (z 2 — z 2 ) 2h 
Taking a conformal transformation (12.481) . for t < l 

(O a (wi,wi)O a (w[,w 1 )O a (w 2 ,w 2 )O a (w 2 ,w 2 )) | 2 —sheets 

1 1 

{w\ — w'-^yb-iuil — W^yb (w 2 — w 2 ) 2 h (w 2 — w 2 ) 2h 
and for l < t < l + L 

(O a (w 1 ,W 1 )Oa(w[,w[)O a (w 2 ,W 2 )Oa(w 2 ,W 2 )) | 2 —sheets 

1 1 1 

da (w\ — w' 2 ) 2h (wi — W^yb (w 2 — w' l ) 2 h (w 2 — W 2 ) 2h 
where we still only keep the leading divergent term for small e. Consequently for l < t < l + L 


+ less divergent term, (2.59) 


+ less divergent term, (2.60) 


dwid'd W2 d / (Oa(w 1 ,wi)Oa(w l ,w 1 )O a (w 2 ,w 2 )O a (w 2 ,w 2 )) 


(2h) 2 (2h + l) 2 


1 


d a 

(2h) 2 (2h + l) 2 


da 


(lU\ — W 2 y h + 2 (wi — w[) 2h (w 2 — W l ) 2h + 2 (w 2 — W 2 ) 2h 

1 

+ lower divergent term, 


(2e) 8/l+4 


+ lower divergent term 


(2.61) 
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which leads to AS 2 = log d a for l < t < l + L. 

Next we consider the descendent operators with the form 

V(w,w) = L^O(w,w), 


(2.62) 


where L ( ^ is a combination of holomorphic generators such that V(w,w) is a quasi-primary 
operatoio. We also assume the operator L has a fixed conformal dimension 


[Lq, )] = mL ( \ 


(2.63) 


Because the final result only depend on the most singular term in the two-point function 
and the 2n-point function, there is nothing change by adding in the operators with smaller 
conformal dimensions. 

The conformal transformation for the descendant operators are different from the one for 
the primary operator. Under a conformal transformation 


L-r 


1 q 

f dw 

2 vr % J 

[w — tt>l ) m_1 

1 q 

( dw 

2m J 

[w — W l) m 1 

1 „ 

f dw 


2 iri 

L_ r 


771—1 


(w — w 1 ) 

(I s -)" + 


£ 


, dz 

'dw 

L_ 


12 

ru '‘ (l i ) 2 + ^u.»}) 


(z — z\) r+2 y dw' 


12 


(2.64) 


-m \z\ vo 

aw 1 

Here the ellipsis denote the terms with lower conformal dimensions, which leads to less di¬ 


vergent terms in the limit e —>■ ojflunder the conformal transformation (12.481) . the four-point 
function transforms as 


(V (wi,wi)V (w^wjV (w 2 ,w 2 )V (w 2 ,w 2 )) 

. dwi . h+m 1 9w\ ^h+m / 9 w 2 n /j_|_ m / dw 2 dw^ h+m / dw^ .. 

[ dz 2 ) { dz 2 ’ y dz[’ ^dzj 

■(V{zi,zi)V{z 1 ,Zi)V{z 2 ,z 2 )V{z 2 ,z 2 )) -l-less divergent terms. 


,dw 2 sh+m ,dw 2 \h+m 

V ^ t J V n.-/ ) 

OZ2 UZ2 

(2.65) 


Even though the operator V(w, w) is not a primary operator, the coefficient for the leading term 
is the same as the one for the primary operator. The terms with lower conformal dimensions 
are less divergent in the small e limit, so do not contribute to the final result. 

Because we assume that the operator V is a quasi-primary operator, it transforms homo¬ 
geneously under a linear conformal transformation. Consider the conformal transformation 


1 


u = — 
z 


( 2 . 66 ) 


3 This condition can be relaxed from the result in the next subsection, but for simplicity we still assume this 
condition in this subsection. 

4 The most divergent term in the limit t —> 0 comes from the OPE of two operators to the identity operator. 
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we have 


Furthermore 


lim z 2{ ~ h+m) L^O(z) 

z— 


= lim {^ h+m h 2{h+m) L^O{u) 
2 ->oo dz 

= (-1 ) h+m L^O{u) | u=0 • 


L—n |it=0 — 


— <f T(u) 
m J u n 1 

T(z)z A 


2 iri 

1 


dz 


2 ni [ z 71+3 


— Lr 


(2.67) 


( 2 . 68 ) 


and 

lim z 2(h+m \0 | L^O(z) = (O | L^-l)™ (2.69) 

Z —>-00 

then the two-point function is 

(L<->0( 2 ), £<->o(*')> = (z _ z y iMm y (2.70) 

where Co = (—1 ) m (h | \ h). 

Now let us study the divergence in the correlation function of four descendants in the z 
coordinate. We set l < t < l + L then {z\ . z 2 ) (z2, z\ ) (zi,^) (z2,z 2 ) are close to each other. 
The four-point correlation function of quasi-primary operators can be transformed into 


(V(z 1 ,z 1 )V(z , 1 ,z 1 )V(z2,z 2 )V(z' 2 ,z' 2 )) 

= V(0(z lt z 1 )0(z' 1 ,z 1 )0(z 2 , z 2 )0(z 2 , z 2 )) 

= 'D'5^ j c m (0(zi)0(z 1 ) \ m 0 (z 2 )0(z 2 ))(0(zi)0(z 1 ) \ m 0{z 2 )0{z' 2 )) 

m 

= 'D'Y2 c m,n(0(zi)0(z 2 ) I m 0(z2)0(z 1 ))(0(z 1 )0(z' 1 ) \ n 0(z 2 )0(z 2 )) 

m,n 

= ^c mj „(L ( - ) 0(zi)L ( - ) 0(4) U L^0(z 2 )L^0(z 1 ))(0(z 1 )0(z' 1 ) \ n 0(z 2 )0(z 2 )), 

m,n 

(2.71) 


where {0(zi)0(z' 1 ) \ m 0(z 2 )0(z 2 )) denote the conformal block expansion with the Virasoro 
module [m] as the propagator. The first equation transforms the correlation function of four 
descendants into the differential on the correlation function of corresponding primaries. In the 
second equation, we expand the partition function by the conformal blocks and c m denote the 
OPE coefficient. In the third equation, we expand the holomorphic part in terms of another 
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channel El- In the fourth equation, we change the differential back into the Virasoro operators 
acting on the primaries in the correlation function. The last equation is based on the fact that 
the Ward identity is satisfied for the conformal blocks. 

The most divergent term only comes from the one with m = n = 0 and Co,o = HU. 
Actually even in the vacuum block, only the identity operator gives the most divergent term, 

{V(z 1 ,z 1 )V{z' 1 ,z 1 )V(z 2 ,z 2 )V(z 2 ,z 2 )) 

= }{LH 0 (, 1 ) L (A0(4)}(LH0( Z2 )L(-)0(4))(0(zd)0(z-;))(0(2 2 )0(4)) 

da 

+less divergent terms. (2.72) 

Changing back into the rc-coordinate and keeping the most divergent term, we find 

(V (wi,wi)V (w 1 ,w' 1 )V (w 2 ,w 2 )V {w' 2 ,w 2 )) 

= 0 (w' 2 )){L^ ] 0 (w 2 )L { ^ 0 (w[))( 0 (wi) 0 (w[))( 0 (w 2 ) 0 (w' 2 )) 

da 

+less divergent terms. (2.73) 

Therefore, for a quasi-primary operator we still have AS 2 = log d a . 

In short, for two kinds of descendent operators: dO and L^O, the difference between 
the second Renyi entropy of the excited states and the ground state is universal, equals to 
the logarithmic of the quantum dimension. This is the same as the case of inserting a local 
primary operator. 


2.5 n-th RE for generic descendent states 

With the above study, we are ready to discuss the effect of most general descendent operators. 
A descendent operator in a module generated from a primary operator O a may take the 
following generic form 


V= ]T d mJ . r ,Ud m L(-^)(d r L(-^)O a (w,w) 

m,j,r,k 


(2.74) 


where < k )o a (w, w) is a quasi-primary operator, and L ( d) (£( > fc )) i s a combination of 

holomorphic (anti-holomorphic) Virasoro algebra with fixed conformal dimension [Lq, Li - ’- 7 )] = 
Keeping the most divergent term in the correlation function, we see that only the 
terms with m + r + pj + pk + 2 h a = A is non-zero. For general n-sheeted surface under the 
conformal transformation (12.4811 


it — l — ie) « 27rw,-_]\ I 

— - () n '' 


(/ H - L — t -K ic) 1 


(2.75) 
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(n) 

Figure 2: The fusion transformations to obtain A S y A 


/ — l ~\~ It) n 27 ri / ■ j ni 

Z- = - r e n O ) n 

(l + L — t — it) 2 


(2.76) 


we see that Zj is close to z'^ +1 as in HU. 

For the n-th Renyi entropy, we need to compute the 2n-point function. Similar to the pre¬ 
vious calculation we need to take a conformal transformation to z coordinate and take proper 
channel to expand the 2n-point function into the holomorphic and the anti-holomorphic part, 
as graphically shown in Fig. [^j]. In each channel only the identity operator contribute to the 
final result, so the 2n-point function breaks up into n two-point functions for the holomorphic 
part (and n for the anti-holomorphic part). We can take a conformal transformation back to 
the w coordinate, with the leading divergent term being transformed homogenously. Actually 
this recipe only needs us to know which two arguments are close to each other and the coef¬ 
ficients for the most divergent term. With this message we can directly do the calculation in 
the w coordinate and multiply the coefficients from the OPE and the channel changing 


(V(w\, w{)V(w 1 ,w 1 )...V{w n ,w n )V(w n , w n )) | n—sheet 




m 1 ,j 1 ;r 1 ,k 1 ' 


e^ ei dm n ,j n -, 


in ■T'ri. ■ h 


rt* 

U r ./ / , / 

‘ ™nJni r n/K 


I n-sheet 


= d~^ E E d 




ct 


' ' k 1 ■■■ E E 


m n ,jn,rn,k n m' n ,j' n y n ,k' n 

(d^L(-’^o a (wyd m h { -^o a (w 2 )) • • • (d^L^-^oywyd^L^Oaiw'y,) 
w 2 w i 

(92 1 Z ( -’ fc Z0 Q (h;i)9:U ( -’ fc i ) 0a(%))...(^Z(-’ fc ™)0 Q (u; n )5;AZ(-’ fc ^0 a (<)) 

w n 

+less divergent terms. 


dm n ,j n ,r n ,k n d r ■/ i ,i 
"°n iJni' nrn 


(2.77) 


Each two-point function in the above relation can be computed directly. For example, for the 
5 We have extracted this figure from im- 
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holomorphic part, it is 


(d^L^^o a ( Wi )d m ; +1 L^'^o a (w' i+ 1 )> 


u i +1 


/ (/i | | hU, ./ (-l)»i 

* fl m *+l •?*>•?»+1 

— U Wi U ’ 


u i +1 


K - u;' 


(h I I/i)<S. ./ (-l)Wi +mi 

(u;* - tw- + i ) 2(ft+w ‘ )+m<+m ‘+ 1 

where we have used (12.701) in the first equation. Introducing the matrices 


(2.78) 


D _ J* 


M { m, j} ,{r,k} = (h | L(-’^t L (-d) | h)6 jtk i r - m , 
and defining the density matrix 

p = BMB^M\ 

then the two-point and 2n-point function can be written as 

(V(w 1 ,w 1 )V(w l ,w 1 )) = ^ 3 A » 

(V(wi,W 1 )V(w[,w[)..V(w n ,W n )V(w n ,W n )) | n-sheet = 
With the normalized density matrix 


trp n 


(2e) 2nA ‘ 


Po = 


P 

trp ’ 


we find 


Z< 


= d-^-'hrp 5 , 


and therefore 


(2.79) 

(2.80) 

(2.81) 

(2.82) 

(2.83) 

(2.84) 


A S n = log d a -5— log trpl = A S" Y -log trp%, (2.85) 

n — 1 n — 1 

and 

ASee = log da - trpo log po = A,S™ ry - trp 0 log p 0 (2.86) 

where 

A S" y = A5gg ary = log d a (2.87) 

is the quantum entanglement of the primary operator. 

It is obvious that there are two kinds of contribution to the entropy of a generic descendent 
operator. The first one takes a universal form, depending on the quantum dimension of the 
corresponding primary operator. Even though the theory could be different such that the OPE 
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coefficients and the channel changing coefficients are different, the 2 ra-point function always 
takes the form (12.771) with a different coefficient, so the relation between the entropies of the 
descendant operator and the corresponding primary operator does not change. 

The extra contribution to the entropy is remarkable. Formally it takes the form —trpo log po- 
It is zero only when the matrix d m j- r j c is of rank one. By matrix product, the matrix po is 
also of rank one. under a conjugate transformation, po can be diagnosed as ( 1 , 0 , 0 ...), which 
has no more correction to the entanglement entropy and Renyi entropy. In this case the en¬ 
tanglement and Renyi entropy are equal to the log of quantum dimension. The two examples 
in the previous section dO(w, w) and L^~’0{w, w ) belong to this class. A general form of this 
kind operator is like 

V(w,w) = jC^£^0(w,w), (2.88) 

where C and C are composed of holomorphic and anti-holomorphic Virasoro generators re¬ 
spectively. 

For the cases when the rank of matrix d m j- r k is more than one, the entanglement entropy 
and Renyi entropy have extra corrections. The simplest example is for di,o ; o,o = d o, 0 ;i,o = 
1 with other coefficients being zero. This gives the operator V = (d + d)0(w,w). The 
entanglement entropy has extra log 2 correction. 

The extra increase for the entanglement entropy is easy to understand in a free theory. For 
the free theory, the entanglement entropy for local operator can be understood by the quasi¬ 
particle. The increase of the entanglement entropy of the local operator from the vacuum is 
equal to the entanglement of the EPR pair [[?[[. The entanglement entropy for the EPR pair 
only depends on the number of the pair and their relative normalizations. For the operator 
(|2.88[> . we change all of the left- and right-moving particles simultaneously. Because all of 
the EPR pairs have the same conformal dimension their normalizations change in the same 
way, which do not change the entanglement entropy. However for a more generic descendant 
operators which can not decomposed as (|2.88|) . it change the relative normalization of the EPR 
pair and even increase the number of EPR pairs. This is the origin of the extra increase in 
the entropies for the descendants. 

From the derivation, we see that only the leading divergent term in the OPE appears in 
the final result. The only thing we should consider is which operators are close to each other 
under the analytically extension 

w —>• —l T (t — ie), w —y —l — (t — ie). (2.89) 

The crucial point here is that the holomorphic and anti-holomorphic parts have different limits. 
For example for the four-point function z\)0-2{z2, 22 ) 03 ( 23 , ^ 3 ) 04 ( 24 , 24 )), the e —>• 0 
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limit set ( 21 , 22 ) close to each other while ( 21 , 23 ) close to each other. The OPE cannot be 
used directly. That’s the reason why we need to transform into the z coordinate and use the 
relation between the conformal blocks in different channel. 

2.6 Comments in BCFT 

In this subsection, we would like to introduce a boundary (or defect) at x = 0 as shown in 
Fig. 1. This boundary preserves the conformal symmetry. There are two kinds of boundary 
preserving conformal symmetry, and they correspond to the Neumann boundary condition 
and Dirichlet boundary condition normally. The global property of the CFT with a boundary 
has been discussed in HZI M M- As we know, the correlation function in a BCFT are much 
different from the one in CFT without boundary. In terms of (12.41) . the Renyi entropy may be 
sensitive to the boundary. In m, the Renyi entropy of the primary states with a boundary 
has been studied and it was shown that the maximal value of Renyi entropy does not change. 
The boundary effects just change the time evolution of Renyi entropy. It is also interesting to 
check what will happen to the Renyi entropy for the local descendent states in 2D CFT with a 
boundary. For simplicity, we only discuss the descendent operator whose entropy is the same 
as the one of corresponding primary operator. 

As shown in mi na, the Renyi entropy highly depends on the conformal blocks of the 
theory for general rational CFTs in 2D. In terms of |23j [2l], the n-point correlation functions 
in 2D CFTs with a boundary are related to the holomorphic part of the conformal blocks of 
the 2ra-point correlation functions on the 2D full complex plane. There is a systematical way 
called the image method to translate n-point correlation functions in 2D BCFTs to 2n-point 
correlation functions in 2D CFT. The more precise relation is that an n-point function in 
the upper half plane(UHP), which is a function of the coordinates ( 21 ,, z n \ 21 ,..., z n ), behaves 
under conformal transformations in the similar way as the holomorphic sector of an 2n-point 
function in the full plane which depends on ( 21 ,..., 2 n ; 2 *,..., 2 *), analytically continued to 
z* = Zj. In El, the time evolution of Renyi entropy highly depends on the holomorphic part 
of conformal block. In 2D CFTs with a boundary, the boundary indeed changes the propagator 
of the primary field but does not change the fusion constants in the bulk. In this sense, one 
can expect that the time evolution of Renyi entropy of the descendent states in CFTs with a 
boundary is almost the same as the one in CFTs on the full complex plane. From the studies in 
the previous subsections, the Renyi entropy of the descendent states could be the same as the 
one of corresponding primary states. One can check that the boundary does not change the 
maximal value of the Renyi entropy of the descendent states. In terms of the quasi-particles 
picture given in m, the boundary just changes the time evolution of Renyi entropy of the 
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descendent states, just as the primary states. 


3 Renyi Entropy in Deformed CFT 


In the 2D rational CFT, we find that the Renyi entropy of local descendent operators could 
coincide with the logarithmic of the quantum dimension of the corresponding primary operator. 
Now we would like to study these operators in the CFTs with additional deformations or 
interactions. In the deformed CFTs, there is no conformal invariance. Nevertheless the effect of 
deformations on the Renyi entropy of local excited states can then be studied within conformal 
perturbation theory [25], if the deformation is weak. Here we only consider the theory which 
is perturbed by the local interaction ^(z). Namely, we consider the theory which is perturbed 
by an operator with the conformal dimension A > 2. The action is now 

I = I 0 + SI = I 0 + X [ $(z)d 2 z, (3.1) 

Jr .1 


where Iq is the original CFT action, A = g/ A 2 ~ a with g being dimensionless and R\ denotes 
the complex plane. Consequently the total Hamiltonian H of the theory is H = Hq+SH, where 
Hq is the Hamiltonian of original CFT. The system will be evolving under the Hamiltonian 
H in the real time approach. As usual, we consider an excited state by acting a primary or 
descendent operator O on the vacuum |0), which is define^ by H |0> = 0. 

We consider the parameter g <C 1 in m ■ The two-point correlation function of the 
primary or the descendent operator is 


(O t (wi,wi)0(a;2,u;2))si = (e $( ~ )cZ w 2 ))f?i, (3.2) 


where Si and R\ denote the complex plane in deformed CFT and original CFT respectively. 
We expand the deformation with respect to the powers of A 

)Ki = g iAA f ... f (3.3) 

N=0 jRl Jri 

The two-point correlation function ([3.211 contains infinite towers of contribution from interac¬ 
tion between the operators O and in CFT. Similarly, the 2n-point correlation function on 
S n can be expanded 

(O t (wi,U;i)0(w2,a)2)...0(w2n,W2n))s n 
= <e -A/^* 0 t (wi ,^ i)0(w2 ,^ ) ... 0(w2n ,^ 2n) > Jin 

= ^ / ••• / (®(zi)...$(z N )0^(uj 1 ,U 1 )0(u2,02)...0(uJ2n,U2n))R n d 2 Z 1 ...d 2 Z N , 

]\[ = 0 ’ ** Rn J Rn 

6 We set the vacuum energy to be zero for convenience. 
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where E n is n copies of Si. To calculate the correlation function on R n , one can employ the 
following conformal transformation 

/ m \ V n c n 

and w = w ® = l eri- (3 - 4) 

Then R n is mapped to the complex plane R\. The points wi, ...,u)2n is mapped to z i, Z2 n 
respectively by (|3.4I) . The 2n-point correlation function can be expressed as 

(O f (wi, wi) O (w 2 , w 2 )... O (w 2 n, w 2n )) s n 

+°° /_\\JV r r 

N =0 7V Jr i jRl 

O t (Cl ! ?l)0(6 5 6)---C ) (6n,6n)) J Rid 2 |l...(i 2 |Af, (3.5) 

where C n is the Jacobian from 0(cui)...0(cu n ). We know that when L + l < t or t < L, 
~ O(e) (i = 1,2 n — 1). Then the leading order term in the 2n-point correlation 
function is 


(O t (wi,o;i)0(a;2,W2)...0(w2n,W2n))s T1 

/_ \\N r r 

- Cn N\ / / ( UJ, (£i)...uj , (£ N )) 2 - A (®(£i)...®(£ N )}R 1 d 2 Z, 1 ...d 2 £ N 

^~ 0 2V! JRi 

n— 1 

JJ ( ot (^2*+i, 6i+i)0(6t+2,6 i+2))iii- (3.6) 

i=l 

In (13.61) . the contribution from 0(zi ) being contracted with 3>(£i) is subleading in the limit 
e —>• 0. There are some subtle issues which have been come across in perturbative CFT. It 
is well known that such integrals are potentially ambiguous due to the singularities from the 
contact terms when 0(zi ) contracts with $(£j) [26| [27] . After proper regularization which 
highly depends on the specific choice of O and T, the contraction between O and make a 
finite contribution. That means the dominant contribution in the limit e —> 0 comes only from 
the contraction between O’s as shown in the last line of (13.61) . In order to present the crucial 
point clearly, we take the second Renyi entropy as an example. The second Renyi entropy 
(12.41) can be expanded by the powers of A as following 


= log 


(oooo)r 2 2 X (oooo ) R2 


<ch( 0 OO) Rl (u/( 0 ) 2 - A d£ 


(OO) 2 


+ 


ao 2 

(OO) 2 


Ri 


(00) Rl 

m)oooo) Rl + 


Ri 


(3.7) 


where the ... stands for the higher order terms in g. In order to make the perturbation well 
defined, we need a regularization to deal with the divergent terms in the integration. Espe¬ 
cially for f R (uj'(^(^) 0(wi)0(w2))R 1 d^, we should introduce two contract terms like 
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f (bo 1 (£)) 2 ~ A (&(i)0(wi)0(w2))R 1 (5(£—wi)+5(£,—W2))d£ : to subtract the divergence if <h is close 
to O. For the second integration in (13.711 , we should do a similar regularization to deal with the 
divergence. All the divergent terms can be absorbed into A to make the perturbation well de¬ 
fined. Finally, {^0{z\)0{z2)0{z^)0{zi)) r 2 can be expressed by {^){0{z\)0{z2)0{zz)0{z^)) r 1 
in the early time limit or the late time limit (12.111) . Repeating the analysis, we can calculate 
the Ato higher order of A formally. For example, 



+ 


log 


(QOOO)r 2 

(ooy Rl 


A C 2 (0000)r 2 


(< 00 ) 2 


Ri 


2A C 2 (0000)r 2 f 

d£(co'(o) 2 - A m)) Rl +..:. 


(3.8) 


In the relation (13.81) . = 1/do which corresponds to the contribution from the local 

\ uu )r 1 

excitation during the time l < t < L + l. 

If is a primary field, (d> )r 1 = 0. We consider the second order of A. The two-point 
function is 


(0 t (wi,wi)0(u;2,u;2)si = wi)0(w 2 , oj 2 )r 1 (1 + ^A 2 f {®(z 1 )$(z2))R 1 d 2 z 1 d 2 Z2). (3.9) 

z JRi 


Similarly, the 2n-point function 




n_1 - 1 r 

Cn JJ(0 t (^2i+l,?2*+l)0(^2j+2,6i+2))(l + / {${zi)<&(z2)) R n d 2 Z\d 2 Z 2 ), 

i=l d Rn 


during the time L + /<iorO<t<Z. Then 


_ 1 + I A 2 / {®(zi)®(z 2 ))R n d 2 zid 2 Z2 

CRT -r "6 (i + i A 2 f($( Zl )<p( Z2 )) Ri d2 Zl d 2 Z2y 

(n) 


A S {n) = AS^lr + Z —— log 
2 ( 1 - 


- as cft + on \ n ) a2 ( J {®(zi) < &(z 2 ))R n d 2 zid 2 Z 2 -n J{^{zi)^(z 2 ))R 1 d 2 z l d 2 Z 2 j, 

(3.10) 


where ASq' rt denotes the result in the CFT. The second part in (13.101) is just the correction 
to the Renyi entropy from the deformation. It is called global contribution in the paper [25l| . 

(n) 

During the time l < t < L + l the only change is in A S y CRT , which contributes the log do- 


4 Conclusion and Discussion 


In this paper, we have studied the Renyi entropy of local descendent operators in 2D CFT, 
extending the previous studies in [8j [Tlj JjjJ. In [8] (FLJ [15], it has been found that the Renyi 
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entropy of a primary state is equal to the logarithmic of quantum dimension of the primary 
operator. It is a natural question to consider the quantum entanglement of the descendent 
states. Firstly, we showed that for the specific operator L_iO with O being primary, its Renyi 
entropy is still the logarithmic of quantum dimension of the primary operator O. Secondly, 
we discuss the descendent state L^'O |0) generated by a quasi-primary operator. For such 
quasi-primary states, we showed that their quantum entanglements are the same as their 
primaries. Despite the fact that the operators look quite complicated and their conformal 
transformations are involved, the leading divergent terms in the early time and late time limit 
are simple, behaving as the one for primary operators. As a result, the quantum entanglement 
of the quasi-primary operators are the same as the primaries. Moreover we discussed the most 
generic descendent operators of the form 

D= £ (4.1) 

m,j,r,k 

Out of surprise, we found that the Renyi entropy of such operator is generally different from 
the one of the primary O a (w, w). Only when the rank of the matrix d m j- r ^ is one, the entropies 
are the same as the ones of the primary. A typical example of such operator is 

V W {w,w) = (4.2) 

Otherwise there is extra contribution. A typical example with extra contribution is the oper¬ 
ator of the form 

V^(w, w) = (L~ + L~)0(w, w). (4.3) 

To clarify the entropy difference between two kinds of operators V W (w, w) and V^ (w, w), 
it would be illuminating to consider the free scalar field theory (j> in 2D. Consider the primary 
operators 0 \ and 0 2 , with 

Oi = e*, 0 2 = -i=(e^ + e-^), (4.4) 

v2 

where cf) = (f>(z) + (j)(z). Following [8] [111 the excited state 0\ |0) is regarded as the product 
state |0) L <8> e l ^ z ^ jO)^ in the chiral and anti-chiral sectors. It is not an entangled state, 
so we get a vanishing entanglement entropy. On the other hand, the operator 0 2 creates the 
maximally entangled state, ■^(e*'^ |0) L (8) |0)jj + e -i< ^ |0) L <8 e~ 1 ^^ |0) R ). The Renyi 

entropy is log 2 when one of the sector spreads into the region of the subsystem. The first 
descendent operators 0[ = L_iOi and O2 = L_i0 2 are 

Oj" 1 ) = id(j)(z)e i< ^, 0 { 2 ~ 1] = J-dfaXe** + e~^). (4.5) 

V2 
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The state 0[ 1) |0) = (id(f>(z)e l ^ z ' > ) |0) L <8 e 1 ^ \0) R is still a product state, has vanish¬ 
ing the Renyi entropy. Similarly, the state o\ ^ |0) = -^=[( idcj)(z)e 1 ^ |0) L ( 8 ) (0)^ — 

(idc/)(z)e~ l ^ z ^) |0) L (8 |0) fi ] is still a maximally entangled state. Therefore, the Renyi 

entropy should be same as the original state. We could do the operation on the chiral sector 
by a ( i,ki = 1,2,...) any times. The result is just a new EPR state. The operation by 

L_j does not change the entanglement property, but just change the state in the chiral sector. 
This is the basic reason that the local excitation by the descendent operators works similarly 
as the primary operator in this class. The same argument could also be used for the operations 
UiL% (i,ki = 1 , 2 ,...). 

Now Let us see what will happen when we apply the operation ^Z{ ni }{ nj } (IT iL-m + 
u, T- n .), e.g., L_i + Z_ 1 . The operators Oi and 0 2 become^ 

Oj -1 ’ -1 ) = (L_ 1 + Z_i)Oi = (id<j)(z) + id^(z))e i ^ +i ^ z \ 
o£ _1 ’ _I) = (L_ 1 + Z_i)0 2 = (id</>(z) + id^{z)){e^ - e~^). (4.6) 

For the operator Oj , the corresponding state is 

oi- 1 ’- i )| 0) = | 2 > L ®| 1)/? + | 1 > L ®| 2 > iJ , 

where we define 

|1)l = |0> t , |2 ) l = |0)i , |l) s = |0) s , |2 ) h = |0) s . (4.7) 

With some normalization this is just the EPR state, the entanglement entropy of which is log 2 
when one of the sector state spreads into the subsystem A. In our calculation this contribution 
comes from the second term in (12.861) . For the operator Og 11 the state is 

oP'-I> |0) = (|2) I |l) s + |1) L |2) a ) - (|4 ) l |3) fl + |3 ) l |4) s ), (4.8) 

where we define 

|3) L = e~ l *U\0) L , |4) i = «^)e- i «W| 0 > 1 , |3) a = e-‘« i )| 0 > a , |4) a = « 0 ( 2 ) e -*»W | 0 ) R . 

We could see that after normalization the relation (14. 8p is just the direct sum of two EPR 
states. The entanglement entropy is 2 log 2 when one of the sector state spreads into the 
subsystem A. This is also consistent with our calculation (|2.86l) . 

Indeed for the interacting theory it is not easy to find such a clear explanation. Actually we 
still need better understanding on the fact that the seemingly decoupled chiral and anti-chiral 
sector are entangled with each other, even for the primary excitation state, see [53] . 

'We ignore the constant -^= for O 2 below. 
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In this paper, we also introduced a boundary which do not break conformal symmetry 
as in [T5] . For the descendent state of type dm we estimated the maximal value of Renyi 
entropy of local descendent states, which is the same as the one in the theory without such 
kind of boundary, although the boundary change the time evolution behavior of Renyi entropy. 
Finally, we deformed the original CFT by additional operator and discussed the Renyi entropy 
of generic local excited states. If the deformation is small, we can treat it as a perturbation of 
CFT. By this way, we can estimate that this deformation make a global contribution to the 
Renyi entropy up to g 2 , if the deformation is a primary field. Though there are divergences in 
the global contribution (253, such divergences can be regulated appropriately. Actually there 
are various definite extensive examples [28] [29] [30] m m to show that global contribution 
is finite. 
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Appendix A Conformal Transformation for Descendent Oper¬ 
ator 

In this section, we show how to take the conformal transformation for a decedent operator. 
For a general primary operator, under a conformal transformation it acts as 

(A.l) 
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First consider the conformal transformation of L_i^(z) 


L-iipiz) 
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2ni 

1 

2ni 

1 

2n i 

1 

2n i 
d 2 i 

W 2 


£ dz z' n+1 T{z')y(z) 


, / rm . t.,dw c , / /,. .,dw. h/ dw.z 

^ )f ar + ist" ■*}]«")< & )( w) 


/ .chn 2 r hipiw) 


dz ( a ' ) [ ( ' \ 2 1 ' 

oh; (re — w) 2 w — w 
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d w y(w) , ,dw. h ,dw. h 


+ normal] (—) (—) 


,dw . 


(^^h+(—)^h)(— r(—) 

dz 


dw h( ^h 

dz’ [ dz 


(A.2) 


This is really the conformal transformation for the operator d'i^(z), which support our calcu¬ 
lation. 

For more generic descendent operator 


L- n ^{z) = j dz (z - z ) n+1 T(z)ip{z) 


1 x , 4 . 1 ,\,dw , 9 c r / . ,dw . h ,dw .t. 

= 2S h l r <“ ><a7> + lA" 


f, > ( > s_ n+ i u dw 2 hip{w) d W i’(w) 

- (pdz{z -z) + (7-7-^2 + —-+ = T{w )ip(w) :) 

1 J dz (w — w) 2 w — w 
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c , / . ...dw.u.dw.z 


(A.3) 


Calculating the contour integral, the conformal transformation of specific descendent operator 
is 
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d 2 i 


d 3 i 


L -^z) = 

dz 


) 2 + i 


dw, h dw - h 




3d 2 w „ ,, ..dw.u.dw.z ,dw. 0 ..dw. h ,dw.z 

+ 5 & )‘( + ( & ) & ) ( fe) 


c , 1 ,dw. 2 




(A-4) 


Actually, we can get a more general conformal transformation for the operators, which is 
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an expansion for the previous result 
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